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Abstract 

The Aero Propulsion Laboratory is currently developing 
a computer-aided design program for high power airborne 
systems. An important part of this design program will be 
the feasibility study which was to be based on summary algor¬ 
ithms. These algorithms were to relate the weight and volume 
of each system component to the contributing operating para¬ 
meters . 

This study first centers on the requirements for an 
interpolation scheme to form the summary algorithms. The scheme 
will need to work in at least four dimensions and produce accu¬ 
rate results over wide ranges. The advantages of the two 
possible interpolation approaches (algorithm development and 
direct interpolation) are described. The results of this com¬ 
parison show clear advantages in the direct interpolation 
approach using stored data. 

The remainder of this study is a description and evalu¬ 
ation of the subroutine which was developed, INTERP. It has 
the flexibility to interpolate a data array with two or more 
independent variables and will output the values of any number 
of dependent variables. The routine also compensates for mis¬ 
sing values in the known data array and issues an error message 
to the user when a test point beyond the data range is input. 
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INTERPOLATING SUBROUTINE FOR 


HIGH POWER SYSTEM DESIGN 

I. Introduction 


Background 

The Aero Propulsion Laboratory has responsibility for 
the development of multi-megawatt, airborne power systems. 
These high-power systems will be required for directed energy 
weapons being developed at the Air Force Weapons Laboratory. 
The components required in these power systems will include 
turbine, generator, transformer, rectifier, and pulse¬ 
forming network. The three steps in the system development 
are (1) the system feasibility study, (2) the detailed com¬ 
ponent design, and (3) the dynamic system simulation 
(Figure 1). 

The first step, system feasibility study, will use 
summary algorithms which relate the weight and volume of each 
component to the contributing operating parameters. These 
algorithms will then be applied over broad parameter ranges 
to determine optimal combinations of components. Operating 
parameters determined in this step will then be converted to 
component design specifications. Computer programs which can 
produce detailed component designs will then be used in step 
two to specify the detailed design for each component. 

The final step will be the dynamic computer simulation 
of the inter-connected components. This simulation will 
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Fig 1. Flowchart of Computer Aided Design Concept 
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evaluate electrical and thermal performance of the entire 
system. The electrical simulation will permit evaluation of 
voltage and current transients and aid in the refinement of 
the control system. Thermal simulation will provide the data 
needed to specify the cooling system requirements. 

The system study is still in the first phase, the 
development of summary algorithms. A summary algorithm des¬ 
cribes the weight or volume of a component as a function of 
operating parameters. Examples of operating parameters are 
input voltage, frequency, and output power. 

Currently, the Aero Propulsion Laboratory is using 
civilian contractors to develop design programs and point 
designs, which cover the expected range of each of the para¬ 
meters. One contract is with the Raytheon Corporation to 
provide 120 minimum weight and volume designs for rectifiers. 

The contract also includes 110 inverter designs and 65 designs 
for pulse-forming networks. These designs will be delivered 
during 1981 and 1982. 

Problem 

The current problem is the need for a computer pro¬ 
gram which can use the data provided in these designs to form 
the summary algorithms needed for system feasibility analyses. 

The first objective of this study was to determine the 
best format for these algorithms. The format will be deter¬ 
mined by the choice of a multi-dimensional interpolation scheme. 
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Scope 


Considerations in this choice will include (1) compati¬ 
bility of the algorithm format with the next step of the system 
analysis, (2) compatibility of the scheme implementation, and 
(3) complexity, range, and accuracy of the scheme. The actual 
writing of the computer program to form these algorithms will 
be the second part of this study. 

Approach 

The first step in this investigation was a review of 
the High Power Studies (a set of technical reports prepared 
for the Aero Propulsion Laboratory) and direct consultation 
with project engineers in the laboratory to consider the ques¬ 
tion of compatibility with the analysis program. The other 
area of investigation was a literature search and consultation 
with Dr. Lee, Chairman of the AFIT Mathematics Department, to 
compare interpolation schemes. 

Sequence of Presentation 

A more detailed analysis of the requirements which 
affect the choice of the interpolation scheme is presented in 
Chapter II. Chapter III describes the two possible approaches 
to the interpolation problem and the advantages of each. The 
development of an interpolating subroutine based on the direct 
interpolation scheme is also proposed. 

Chapter IV describes the subroutine developed, INTERP. 

In Chapter V the accuracy and effectiveness of INTERP is evalu¬ 
ated using the program CGEN and stored data based on conventional 
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generator designs. Conclusions and recommendations for future 
study are presented in the final chapter. 

There are also two appendices. Appendix A includes a 
"User Guide" to INTERP along with a sample program and a 
listing of the subroutine. Appendix B contains a listing of 
CGEN as well as some of the inputs, outputs, data, and graphs 
used in the evaluation of INTERP. 
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II. Investigation of Requirements 


As stated earlier, the primary application of the sum¬ 
mary algorithms will be in the computer program used to deter¬ 
mine the optimal configuration of system components. This 
optimal design will include some, but not necessarily all, of 
the possible system components indicated below with their 
associated parameters: 

1) Turbine (including gas generator and fuel supply)-- 
speed, power, run-time, number of starts 

2) Generator (three different types will be considered) 
--speed, voltage, power 

3) Transformer -- voltage in, voltage out, power, 
frequency 

4) Rectifiers -- power out, frequency, voltage out 

5) Inverters -- voltage in, voltage out, power out, 
frequency 

6) Pulse-forming networks -- power, pulse width, 
pulse repetition time. 

Figure 2 depicts a possible system design. 

Summary algorithms -have already been developed for 
transformers, using careful analysis of a large number of 
smoothly varying design points. Some algorithms have been 
developed for generators. Due to discontinuities, they are 
valid only over a limited range. Below is an example of one 
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of these algorithms (Ref 2:4). 

P 44Q 

Specific Weight (lbs/kW) * .157 [1.28 - .28 (-gO ’ ] 

. [-.06 + 1.06(^^)'* 6205 ] • [.8567 + .1433(^)J (1) 

where 

P = generator output power (megawatts) 

RPM = rotational speed 

V = generator output voltage (volts, line-neutral) 
Algorithms developed so far for power conditioning components 
are based on preliminary designs and give unacceptable results 
in most cases. 

New design points, now being developed by Raytheon, 
will probably provide the most severe test of the interpolating 
scheme. Large data sets over wide ranges will be developed. 
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Inverter designs based on four independent parameters, with 
some discontinuities, will be included. This will require an 
interpolation scheme, in at least four dimensions, which can 
produce reasonable accuracy over wide ranges. 
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III. Theoretical Approaches 


There are two possible approaches to an interpolating 
problem. The first is to use the available data to generate 
an algorithm, similar to Eq (1), which can then be used inde¬ 
pendently of the data. The alternative approach is to retain 
the data in an array or file. When a value of the function 
is needed, it is necessary to find the nearest known points 
in the array and estimate the value of the function at the 
new point. The next step is to compare these two approaches 
and determine which is best suited to this application. 

Algorithm Approach 

First to be considered is the algorithm approach. This 
approach has been used to date by engineers working on the high 
power systems under consideration. It has the obvious advan¬ 
tage of succinctness. When incorporated into the overall 
design feasibility program, the algorithm approach allows 
faster execution and requires much less core memory. However, 
writing a computer program to reduce a multi-dimensional array 
of data to a compact algorithm proves to be a very complex task. 

One way to design such a program would be to use the 
multi-dimensional spline technique discussed by de Boor in 
Reference 3. At first this approach appears reasonable, 
though quite involved. A problem appears near the end of the 
development where a system of perhaps several thousand 
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simultaneous equations must be solved. Even though the system 
has considerable structure, so that highly efficient solution 
schemes could be used, this is a difficult problem. An alter¬ 
native to the spline approach is to use an application of the 
multi-dimensional Taylor series expansion (Eq 2). 


F(x,, x 7 , 


x„) = F(xi, xi, . 


X A> 


, 3F ( , . 9 F ( 

* Cx l • x l ) * Jx 2 < x 2 • x 2> 

3 F 

.. (x - x^) + higher order terms 

n 


( 2 ) 


where 

(x^, x^, ... x^) = point where the value of the 
function is known 

The key consideration in this approach is the "higher order 
terms." For a nonlinear function there is no limit to the 
number of higher order terms, each of which becomes more com¬ 
plicated. The functions being considered will definitely not 
all be linear. However, a function can sometimes be made 
linear using a change in variables. For instance, if a 
function F varies as x + 5, we could substitute X’ = x +5; 
and F would be linear with respect to x'. This is the approach 
which has been used on several data sets in the past. For data 
which we can examine and manipulate by hand with the help of 
human intuition, we have a reasonable likelihood of discover¬ 
ing variations like ax^ + b, ax c + b, e x + a, or a x + b. 


Physics may also suggest these groups. However, to program 
a computer to uncover all of the possible variations would be 





I 


a gigantic, if not impossible, task. An additional problem 
is that these algorithms are based on the behavior of the 
function near the midpoint of the data sets and often produce 
serious inaccuracies when applied to points near the edges 
of the data array. 

Direct Interpolation 
Appro a~cK 

After considering the complications associated with 
the algorithm approach, a careful look at the direct inter¬ 
polation approach is definitely warranted. An advantage to 
this approach is that output is limited to the range between 
the neighboring data points. This avoids the possibly extreme 
errors of a modeling approach. If more accuracy is required, 
a spline or other more powerful interpolation technique can 
be used. 

As mentioned previously, this approach has the disad¬ 
vantage of requiring more time and memory when being executed 
in the system feasibility program. The objective of the 
feasibility program, however, is not to run in minimum time 
or with minimum memory. Whether it requires ten minutes and 
15K of memory or six hours and 150K, in either case, if the 
program is successful, it will save weeks or months of research 
and greatly reduce the acquisition time of the final system. 

Another disadvantage to be considered is the require¬ 
ment to store and maintain the data sets. This is normally 
done on a tape at the computer center and requires careful 
documentation to ensure that the tapes are retained and that 
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the data on them is readily accessible. The other advantages 
to be gained certainly warrant the small additional effort 
required to document these data sets. 

Conclusion 

An objective of the program developed, using this 
interpolating scheme, will be its applicability to operate on 
data sets still under development. To meet this requirement 
and to ensure reasonable accuracy from the scheme, it is con¬ 
cluded that the direct interpolation approach, using stored 
data, is the best for this application. 

The next step in the study is the development of a 
general subroutine to accommodate at least five variables. 

This subroutine should be general enough to be used on a data 
set from any of the high power system components. Inputs to 
this routine would include the number of independent vari¬ 
ables, the number of data points for each variable, the values 
of the variables at the known data points, an array containing 
the function values at the data points, and the point at which 
the function is to be evaluated. An advantage to this scheme 
is that three values (weight, volume, and efficiency) can be 
stored in the same data array, and all three can be determined 
with one pass through the subroutine. 





IV. Description of SUBROUTINE INTERP 


SUBROUTINE INTERP was developed by the author to meet 
the requirements discussed in the previous chapter. The sub¬ 
routine can accommodate two or more independent variables and 
is not limited to five. It can store any number of dependent 
values at each point in the array. This allows storage and 
interpolation of weight, volume, efficiency, and other data. 

The routine includes an error code output which sig¬ 
nals the user when he has requested data from a point outside 
the stored array bounds. If data are unavailable for some 
points in the array, the user fills these locations with zeros 
and the routine automatically searches until non-zero values 
are found on which to base the interpolated values. Since 
these values normally represent values of weight, volume, or 
efficiency, a value of zero should not otherwise appear. 

In the following section of this chapter is a dis¬ 
cussion of the inputs to the routine and how the data arrays 
are stored. Subsequent sections give an overview of the 
routine and then describe each of the three main sections in 
more detail. A complete listing of INTERP, a user guide, 
and a sample program are included in Appendix A. 

Arguments 

The first input to the routine, NDIM, is an integer 
which specifies the number of independent variables or 

f 
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dimensions to be handled. Allowable values are two or greater. 
The second input, LV, is a vector (length equal to NDIM) which 
contains integers corresponding to the number of data points 
in each dimension. These points do not need to be evenly 
spaced. The order of these points establishes the order in 
which the dimensions will be considered throughout the routine. 

LVT is the next input and is the sum of the elements in 
the length vector, LV. LDI is the number of dependent values 
stored at each point in the array. Any positive integer may 
be used. LDJ is the product of the elements in vector LV. 

LVT, LDI, and LDJ are used along with NDIM to dimension the 
remaining arrays. 

The next input is the vector VAR with length equal to 
LVT. The values of the independent variables at the known 
data points are stored in this vector. The values of the 
first variable are stored first in increasing order, followed 
immediately by the values of the next variable stored the same 
way, and likewise through the independent variables. These 
variables must be stored in the same order as the length values 
stored in LV. The input WVE is the array containing the depen¬ 
dent variables (typically weight, volume, and efficiency). 

The array has dimensions LDI by LDJ. The first index of the 
array specifies which dependent variable to consider for a 
given set of independent variables. The second index walks 
through the possible data locations by first incrementing 
through all values of the first independent variable, then 
through all possible combinations of the first and second 
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independent variables, and so on through all the possible com¬ 
binations of independent variables. Zeros should be stored at 
any indices where the data are not available. The routine will 
check for zero values and search for the next non-zero value. 

The input XV is a vector (length equal to NDIM) which 
contains the coordinates of the test point, the point in the 
data array where the values of the dependent variables are to 
be estimated. The vector LWK (length equal to NDIM) is a work 
space providing temporary storage locations for index values 
within the routine. OUT is the output vector of length LDI 
containing the computed values of the dependent variables. 

The final argument of the subroutine, IERR, is a two- 
digit error message output. The least significant digit will be 
either zero or nine. Zero indicates that a requested point is 
less than the lowest point, and nine indicates a requested point 
greater than the highest data points. The leading digit indi¬ 
cates the dimension in which the out of range condition occurred. 


Mathematical Basis 

The mathematical basis for INTERP is the multi¬ 
dimensional Taylor series expansion introduced in Chapter III 
and repeated below in Eq (2). 


F( x X » x 2> ••• x n ) = p ( x i» x 2’ ••• x n^ 

♦ M, * HL <*2 - 


3F 

+ ‘' 3x X n ‘ x n^ + higher order terms 


( 2 ) 
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The approach used in this routine is a first-order approach 
and neglects the higher order terms. The values of the inde¬ 
pendent variables, x^, ... x n> at the point where the 

function is to be evaluated are labeled XV(1), XV(2) , ... XV(n) 
in the subroutine. 

The subroutine is divided into three main sections as 
shown in Figure 3. The first step in applying Eq (2) is to 
search the data to find the indices of the dependent and inde¬ 
pendent variables at a nearby known data point to be used for 

(xi, xi, ... x'). The next section of the routine checks for 
1 2 n 

missing values (zeroes) in the dependent data locations speci¬ 
fied above and increments or decrements the indices to find 
the nearest, non-zero, dependent variables above and below the 
test point. 

In the final section of the routine the partial differ¬ 
entials are approximated by first difference expressions which 
are evaluated in each dimension and then summed to find the 
desired output. In the remainder of this chapter, these 
three sections are discussed in more detail. 



Fig 3. Flowchart of Main Sections of INTERP 
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Finding Index 

The first section of the routine contains two main 
loops, each with an internal loop as seen in Figure 4. The 
first main loop searches each dimension, or string of indepen¬ 
dent variable values (VAR), to find the index of the lowest 
value which is greater than or equal to XV(I). When this value 
is found, it is stored in LWK(I). If XV(I) is less than the 
lowest value or greater than the highest value in the string, 
an error code is formed and printed, and control returns to 
the main program (note paths A and B in Fig 4). The second 
main loop of this section uses the computed values in LWK and 
the values in the length vector, LV, to compute INDEX. INDEX 
corresponds to the second index of WVE (the known dependent 
variable array) at a point near the desired value specified 
in XV. 

The error indicator, IERR, is initialized to zero at 
the beginning of the routine and is revalued when necessary 
by the output error code as shown in Figure 5. The most signi¬ 
ficant digit of IERR comes from I, the loop index, and indicates 
the dimension being searched when the out-of-range condition 
occurred. The least significant digit is left zero in the case 
of a below range condition (path A), and is changed to nine 
in the case of a value of XV beyond range (path B). 

Compensation for 
Missing Data 

Section 2 of the routine is actually within two do 
loops of section 3 (see Figs 6 and 7). The discussion of both 
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Fig 4. Flowchart of Section 1, Finding Index 
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Fig 5. Flowchart of Output Error Code 

sections will be clearer if section 2 is considered first. 

It must be recalled that TMP2, IT1, and IT2 are already 
assigned in section 3 and that section 2 is included within 
loops with indices K and L. 

The source of the data used as inputs to this routine 
are usually the result of point designs developed by the Aero 
Propulsion Laboratory or by civilian contractors. Examining 
existing data reveals that in many cases data for some of the 
points within the array are not available. In this routine 
the missing values of WVE are filled with zeroes. Since these 
values normally represent values of weight, volume or effici¬ 
ency of a system component, a value of zero should not other¬ 
wise appear. The final section of the program will use the 
values of the dependent variables at points above and below 
the test point in each dimension. This requires that many 
points be checked for zero. 
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Fig 6. Flowchart of Section 2 
Compensation for Missing Data 


















Fig 7. Flowchart of Section 3, Difference Expressions 
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In the first part of this section the index of the 
point above the test point, INDHI, is initialized as index 


i 

! 

i 

i 

i 

i 

] (see Fig 6). The stored value of WVE (K,INDHI) is then tested 

; for zero. If a zero is found, INDHI is incremented and the 

check is made again. If the check fails after the highest 
value in the string has been checked, an error output is gen¬ 
erated (path b) as in section 1, indicating that the point is 
beyond range and in which dimension. Note that the index of 
the independent variable, originally stored in LWK, is stored 

t 

in LHI and is incremented along with INDHI. This value is 
1 required in the final section of the routine. 

( The remainder of this section repeats the process just 

i described for the known point below the test point. In this 

case, the dependent variable index is stored in INDLO and the 
dependent variable index in LLO. Again, a route to the error 
code is provided (path A) when a non-zero value is not found 
in the string. 

Difference Expressions 

The objective of the final section is to determine the 
estimated values of the dependent variables at the test point. 
This is done by approximating Eq (2), using first difference 
equations to replace the partial differential. TMP1 divided 
by VAR(IT2+LHI) - VAR(IT2+LL0) replaces the partial differen¬ 
tial and (x n = x^) is replaced by XV(L) - VAR(IT2+LHI). Com¬ 
pare Figure 7 and Eq (2). The quantity WVE(K,INDHI) in the 
s final statement corresponds to F(xJ, x£, ... x^) in Eq (2). 
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The results of these approximations are discussed in 
the following chapter- It is shown there that the dependent 
variable values are output unchanged when the test point coin¬ 
cides with a known data point, and that the output corresponds 
to a straight-line approximation for points between known data 
points. 
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V. Evaluation of SUBROUTINE INTERP 


The objective of the evaluation was to choose a repre¬ 
sentative data base and develop the program to call INTERP. 

The output of this program is then studied to ensure that the 
subroutine is producing reliable data. The data set selected 
contains the results of a conventional generator design pro¬ 
gram. 

This data set is well suited for a test of the routine 
with three independent variables and three dependent variables. 
The independent variables are RPM, voltage, and power. The 
dependent variables are specific weight, efficiency, and 
volume. This data set has the added advantage that an algorithm 
has already been developed to approximate the variation of spe¬ 
cific weight as a function of the three independent variables. 
This algorithm is the one presented earlier as a representative 
case of three algorithms (Eq 1, page 7). This will provide 
the opportunity to compare the subroutine output with the 
algorithm results. 

Program CGEN 

The program developed to manipulate the conventional 
generator data is named CGEN. A complete listing of the pro¬ 
gram, along with input and output data, is located in Appendix 
B. As seen in the flowchart in Figure 8, the program begins 
by printing the labels for the output data. Both the dependent 

24 


r 







Fig 8. Flowchart of Program CGEN 
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and the independent variable values at the test points are 
included in the printed output. The dimensioning information 
and values of the independent and dependent variables at the 
known points are then read from the data file. As anticipated, 
there were numerous cases where the dependent variable values 
were not available. Zero fill was used at these positions, 
as previously discussed. 

The remainder of the program reads a test point, XV, 
and calls INTERP. If no error is encountered, the output is 
printed. The program then computes the specific weight using 
the algorithm presented earlier (Eq 1) and prints that value 
for comparison. If an error is encountered, the message will 
be printed by the routine and the program goes on to the next 
test point. The program stops when it has used all the test 
points in the input data file. 

Results 

Many test points were run through CGEN to ensure that 
the routine was working properly. Several errors in the 
routine were discovered and corrected before it reached its 
present form. Test points were run at known points to ensure 
the routine agreed there. Points above and below the range 
were used to test the error code function. Points were also 
tried at and near locations with zero fill to ensure section 
2 of the routine performs correctly. 

With all corrections made and all functions of the 
program and routine working correctly, the test points included 
v* 
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Fig 9. Comparison of Algorithm and Subroutine 
Changing Voltage, Volts .65 

in Appendix B were used to produce a set of data for closer 
study. In one group of test points, the voltage was varied 
while holding the RPM at 18,000 and the power level at 8.83 
megawatts. The results are plotted in Figure 9. As can be 
seen in the figure, the output from INTERP agree with the 
given data at the known points and form a straight-line 
approximation between known points. The value from the algor¬ 
ithm is close to the given data at the lowest voltage, an 
error of 6.6 percent. However, at the highest voltage the 
error increases to 40 percent. 
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There are two additional groups of data presented in 
Appendix B. In one group the voltage is held constant at 
577 volts and the power at 4.91 megawatts. The RPM is varied 
from 8,000 to 20,000. The subroutine output again agreed with 
the known data and a straight-line approximation between points. 
The algorithm data agreed much better in this case with the 
error varying between 1.2 and 12 percent. 

In the other portion of the data set, the RPM was held 
at 18,000 and the voltage at 2885 volts, while the power level 
was varied. The subroutine output again agreed well with the 
data as expected; however, the error in the algorithm output 
was nearly constant at 40 percent. 

Plots of all three sections of the data are included 
in Appendix B. Any number of data sets and plots could have 
been listed and plotted, but these three were chosen to show 
the range of possibilities and advantages of the subroutine 
approach over the algorithm approach used previously. 


< 
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VI. Conclusions and Recommendations 


The original objective of this study was the develop¬ 
ment of a computer program to produce the weight and volume 
algorithms required for the computer aided design of high 
power systems (Fig 1, page 2). These algorithms were to be 
based on the data contained in point designs being produced 
by contractors. A careful comparison of the advantages and 
disadvantages of an algorithm program versus a subroutine to 
assess and interpolate the data from these designs was made. 

The result of that comparison was the decision to develop an 
interpolating subroutine. 

Conclusions 

INTERP is the subroutine that was developed. It has 
the flexibility to interpolate a data array with two or more 
independent variables and will output the interpolated values 
of any number of dependent variables. Missing values in the 
stored data are compensated for by using zero fill at these 
data locations. Output values at known points agree exactly 
with the stored values, and the values between points are based 
on a straight-line approximation in each dimension. 

The results of the evaluation made show that the rou¬ 
tine meets all the objectives. The flexibility of the routine 
exceeds the original requirement to handle five independent 
and three dependent variables. Accuracy is as expected, and 
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an error code informs the user when a desired test point is 
outside of the data range. A comparison with a previously 
developed algorithm shows that INTERP can greatly reduce 
errors that sometimes exceeded 40 percent using the algorithm 
approach. 

Recommendations 

One of the next steps in this computer-aided design 
program will be the development of an executive program to 
implement the design tradeoff study in Figure 1. This program 
will iterate through large numbers of possible component con¬ 
figurations to find designs which optimize weight, volume, 
and efficiency. 

Several of the component design programs have already 
been developed and a major task will be the interfacing of 
these programs, available data sets, and the subroutine 
INTERP. Another project that may prove beneficial is the 
extension of INTERP to give it the capability to also extra¬ 
polate beyond the given data range. It is not yet clear 
whether there will be a requirement for this capability. 


V. 
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APPENDIX A 


SUBROUTINE INTERP 









User Guide to INTERP 

This routine is used to interpolate between known 
data points within a data array containing two or more inde¬ 
pendent variables and any number of dependent variables. The 
user must supply the routine with known values of both the 
independent variables (VAR) and the dependent variables (WVE) 
at some known points. Values for the independent variables 
are then input (the test point) and the outputs of the routine 
are the dependent variable values associated with the test 
point. 

If a value of the test point is outside the range of 
the known values, an error message will be output, indicating 
which dimension (or which of the independent variables) was 
out of range, and whether the test value was below or beyond 
the data range. This error code (IERR) is further discussed 
in the argument description below. 

To call INTERP, the following statement is used: 

CALL INTERP (NDIM, LV, LVT, LDI, LDJ, VAR, WVE, XV, 

LWK, OUT, IERR) 

Arguments must be dimensional and valued as described below: 

NDIM -- Number of dimensions or independent variables 
(integer, two or higher) 

LV -- Length vector containing number of known data 
points in each dimension (integer vector, 
length equals NDIM) 
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LVT -- Sum of the elements in the vector LV (integer) 
LDI -- Number of dependent variables (any positive 
integer) 

LDJ -- Product of the elements in the vector LV 
(integer) 

VAR -- Values of the independent variables at the 
known points (vector, length equals LVT) 

WVE -- Values of the dependent variables at the known 
points (array, LDI by LDJ), use zero fill for 
array points with no data 
XV -- Test point, coordinates where the dependent 

variables are to be estimated (vector, length 
equals NDIM) 

LWK -- Working space (vector, length equals NDIM) 

OUT -- Output, containing computed values of the 

dependent variables (vector, length equals LDI) 
IERR -- Error message (output), zero indicates no error 
otherwise, first digit indicates dimension, 
second digit indicates out of range; zero - 
below data range, nine - above data range. 


Sample Program 


Program SAMPLE is listed below as an example to the 
user. The known data is shown in Table I. Notice there are 


two independent variables (VAR1, VAR2) and one dependent vari¬ 
able (WVE). If there was a second dependent variable, its 
values would have been stored in WVE (1,7) through WVE (1,12). 
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C. 


The test points used and the resulting output are shown 
in Table I for comparison with Table I. If one of the dependent 
variable values listed in Table I had not been available, it 
could have been replaced by a zero. INTERP would then have ig¬ 
nored the point and made the interpolation based on the remain¬ 
ing points. 

TABLE I 


Input Data for Program SAMPLE 


vNvari 

VAR2\V 

5.0 

10.0 

20.0 

20.0 

1.0 

2.0 

3.5 

50.0 

3.0 

4.0 

7.6 

WVE 


TABLE II 


Test Points and Resulting Output (SAMPLE) 


VAR1 

VAR 2 

OUT 

6 

25 

1.53 

8 

. 30 

2.27 

12 

40 

3.35 

16 

45 

5.47 

18 

55 

* 


*Error (IERR=29), VAR2 Beyond Range 
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o o o o o o o o o o o 


PROGRAM SAMPLE IMPUT* OUTPUT) 

DI MENS I ON LIVE < 1 - 6 > » VAR < 5 > , XV < £ > * LWK •:£> , OUT < 1 > * LV 


DIMENSION INFOPMRTION 

NDIM=£ 

LV <1)=3 
LV<£>=£ 

LVT =5 
LDI = 1 
LDJ=6 


KNOWN INDEPENDENT VhP I ABLE VALUES 


•‘AR 

' 1 > 

=5 


•'AR 


= 1 

o 

•'AR 

:'3> 

=£ 

o 

•'AR 

-4 :• 


0 

•‘AR 

: 5> 


0 


KNOWN DEPENDENT VARIABLE VALUES 


LIVE 

1 * 1 > 

LIVE 

1*£> 

LIVE 

1 * 3> 

LIVE 

1 * 4) 

LIVE 

1*5) 

LIVE • 

1 * 8 > 


READ TEST POINT 
15 READ** ' XV - I') * 1 = 1 * NDIM> 

IF EOF •: 5LI NPUT> . NE. 0. > STOP 
CALL INTEPP CNDIM* LV* LVT* LDI * LDJ* VAR* LIVE* 
1 XV * LLIK * OUT * IEPP> 

PRINT** "OUT*".OUT <T> 

GO TO 15 
END 


Input Data 


Output Data 


£ £5 
8 3 0 

12 40 
16 45 
18 55 
♦EOP 
♦EOF 


OUT=l.533333333333 
OUT=8.£££688888867 
□UT=3.353333333333 
OUT=5.476888888,88? 
IEPP=£9 

OU T=5. 4 76666666£ 6 7 

STOP 
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Listing of INTERP 



SUBROUTINE INTEPPCNDIM*LV* LVT* LDI*LD J*VAP*WV£»XV* 

1 LUK* OUT * IEPP.J 

DI MENS I ON MVE (LDI*• V:=lP CLVT > »XV CND I M> » 

1 LM!< •: Nil IM) * OUT C L DI * LV CNDI M> 

C SECTION 1. FINDING INDEX 

IEPP=0 
MIN=1 
MAX»LVC1> 

C FIND LWK TO FORM INDEX 

r- 

DO c'O 1 = 1* NDIM 

C CHECK FOP VALUE BELOW RANGE? IF SO* GO TO ERROR CODE 

L 

IF ( XV CI > . L T. VAR CM I ,N> > GO TO 3 0 
LST=LV C I> M 

r- 

C LOOP TO FIND THE MINIMUM STORED VALUE* WHICH IS GREATER 
C THAN OR EQUAL TO THE TEST POINT 

DO 5 K = 1* LST 
J=K-1 

IF CXV C I> . LE. VAR CM IN «-J> > GO TO 10 

C CHECK FOR VALUE ABOVE RANGE* IF SO* GO TO ERROR CODE 

IF CXV CI> .GT. VARCMAXV) SO TO 35 
5 CONTINUE 


C STORE THE INDEX OF THE INDEPENDENT VARIABLE FOUND IN LWKCI) 
10 LWK CD =1 f.J 

C SET MAX AND MIN FOR NEXT DIMENSION 

MAX=MAX*LV Cl «-l> 

30 MIN=MIN+-LV CI> 

GO TO 33 

C OUTPUT ERROR CODE 


30 IERR=I *10 

31 PRINT »*"IEPP="* IEPP 
RETURN 

35 IERR=I *1 0*-9 

GO TO 31 
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o o n o o o o n n o o o o o o o o o n r 


C FIND INDEX OF MVE CORRESPOND I No TO LMX 

39 INDEX= 0 

DO 50 N=£« NDIM 
ITENF=1 
NDE=N-1 

DO 40 M=1.NDE 

4 0 I TEMP=I TEMP *LV <M > 

5 0 I NDEX= (LW!< (N> -1 :< * I TEMP (-1NDEX 

INDEX =INDEX +■ L Mfc < 1 > 

SECTION 3* LOOPS TO FORM OUT<K> 

DO 70 :< = 1 * LDI 
TMP£=0 
I TE'=0 
IT1 = 1 

DO SO L=1jNDIM 

SECTION E> COMPENSATION FOP MISSING DATA 

INITIALIZE I NDH I S-CHECK FOP ZEPO AT MVE C< > INDHI> 

INDHI=INDEX 
LHI=LWfc<L> 

MMAX=LV<L>—LMK<L> 

DO 55 M=1 * MMAX 

IF (MVE (K ' INDHI > . NE. 0) GO TO 56 
IF ZEPO* INCPEMENT INDHI AND P IECK 
LHI=LHIM 

i I NDH I = I NDH I ■»-1T1 

IF END OF DATA STRING IS PEACHED * SEND EPPOP MESSAGE 
I=L 

GO TO 35 

INITIALIZE INDLO SCHECK FOP ZEPO AT MVE Of»INDLO) 

INDLD=INDEX-IT1 
LLO=LMI- <L> -1 
NMAX=LLO 

DO 57 N=1»NMAX 

IF ( <MVE (>: « INDLO) .NE. 0> .'OP. <XV<L> . EG. VAP (IT£*LHI> O 
GO TO 

IF ZEPO DECREMENT INDLO 

LLO=LLO-l 
INDLO=INDLO-IT1 


57 
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I IF beginning OF string is perched, send error message 


I=L 

S3 TO 30 

C form the difference expressions AND sum them to get OUT oo 

SP TMP1 =.MVE • INI'HI > -WVE <K. IND1_3> 

TNP2=TNP2+ < <XV <’L> -V*P<IT£*LM» *TMP1/ 

1 ( VHP < I T£ 4-LHI > -VHP < I T£ +LL3> > > 

r RESET TEMPORARIES IT1 AND IT£ FOP NEXT DIMENSION 

IT1=IT1*LV<L> 

60 IT£*IT£*LV<L> 

7 0 OUT (K - ) =!>iVE •'!< « INDHI .) *-TMP£ 

RETURN 

END 


,r 































*No Data Available 
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<■» 


j Independent Variables 

Dependent Variables 

RPM 

(* 1000) 

Voltage 

(volts) 

Power 
(MW ) 

Sp. Wt. 
(lb/kW) 

Efficiency 

(%) 

Volume 
(cu. in 
* 10001 

8 

2885 

6.87 

. 215 

95.4 

18.8 

10 

2885 

6.87 

.186 

95.8 

12.2 

12 

2885 

6.87 

.173 

95.8 

9.3 

14 

2885 

6.87 

.154 

• 96.0 

7.6 

16 

2885 

6.87 

.151 

95.6 

6.7 

18 

2885 

6.87 

.142 

95.5 

6.0 

20 

2885 

6.87 

* 

* 

* 

8 

1154 

8.83 

.188 

95.7 

19.1 

10 

1154 

8.83 

.168 

96.0 

12.8 

12 

1154 

8.83 

.153 

95.9 

9.9 

14 

1154 

8.83 

.136 

95.9 

8.1 

16 

1154 

8.83 

.138 

95.4 

7.5 

18 

1154 

8.83 

.131 

95.1 

6.8 

20 

1154 

8.83 

* 

* 

* 

8 

577 

8.83 

* 

* 

* 

10 

577 

8.83 

* 

* 

ft 

12 

577 

§.83 

* 

* 

* 

14 

577 

8.83 

* 

* 

* 

16 

577 

8.83 

* 

* 

ft 

18 

577 

8.83 

* 

* 

* 

20 

577 

8.83 

* 

* 

* 


• *No Data Available 


* 
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o o o o o o o o o o o o o n o o o o o o o o o o 


Listing of Program CGEN 

PROGRAM f sen ■ I NF'UT. OUTPUT:' 

Ii I MENS ION UVE <3* 1 05') ,VAPfl5) , XV <3> *LUK<3> < OUT <3) . LV C3) 

C 

C PRINT COLUMN HEADINGS FOP OUTPUT DATA 
C 

£ PR I NT♦ * "PPM *k V L-N PWR <MW) SF'WT LE/RW EFF V 

iin" 

4 FORMAT <2 CF5. 0*5X> * 4 CFG. 3* 6X) > 

READ INFOPMAYION NEEDED TO DIMENSION ARRAYS 

READ** NDIM * LDI ?LD J* L’v'T 
READ** 'LV'IJ* 1 = 1* NDIM) 

READ THE INDEPENDENT VARIABLE VALUES ‘ 

PE AD* * CVAP c I :■ * I = 1 * LVT> 

READ THE DEPENDENT VARIABLE VALUES 

DO 10 J=1 *LDI 

1 0 READ** < I.IVE ' J* k > * K= 1, 1 05) 

READ THE TEST POINT 

15 READ* * CXV •: I > , I = l, ND IM> 

IF cEOF C5LI NF'IJT :• . NE. 0. > STOP "END OF RIJN " 

CALL INTERP <NDIM*LV*LVT*LDI* LD J* VAR* WVE* 

1 XV* LWK* DUT * IERR) 

IF THERE WAS AN ERROR* GO TO THE NEXT POINT 
IF <IERR.NE.0) GO TO 15 

PRINT THE INDEPENDENT AND DEPENDENT VALUES AT THE TEST PDINT 

PR I NT 4 * •: CXV < I > * I = 1,3> * <OUT C.J) , J= 1,3>> 

COMPUTE THE SPEC. WT. USING THE ALGORITHM* FOP COMPARISON 

:Sl.i.l=. 157* <1.28-. 28* < <XV •: 3 :•5 :> . 449') ) ♦ •: -. 08+1. 06* < CXV C 1 > / 14) 

►»'■-. 6£ 05-* ) * C. 8567+. 143:3* CXV •:£> ♦S-G'PT <3. ) 1 000) > 

PRINT THE SPEC. WT. COMPUTED 

F'R I NT* * " SP WT ■: FROM ALG > = " > SW 

GO TO 15 

END 
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Output from CGEN 


105.15 


8 * 10? 

12.14« 

16. 1 

8. 2 0« 

577. 11 

54. 1 

731.2 

3 03.2 

885. 4 

91.6 

87. 

. 224 

. 194 

16 0 

. 140 

. 1 34 . 

124 

. 122 

.231 

. 135 

164 

145 

. 128 

. 132 

226 

. 195 

. 162 . 

157 

. 146 

. 133 

. 00 i'i 

837 

195 

. 160 

. 149 

137 

. 134 

.245 . 

197 

.17 0 

. 160 

. 149 

147 

142 

. 172 

. 154 

139 

. 137 

.121 . 

120 

. 197 

. 175 

. 154 

140 

134 

. 0 00 

.201 

183 

. 156 

.149 . 

138 

. 130 

. 00 0 

. 000 

177 

16 0 

. 143 

. 135 

0 0 0 

.215 

. 186 . 

173 

. 154 

. 151 

. 142 

0 0 0 

0 0 0 

. 000 

. 000 

0 0 0 

. 000 

. 000 . 

188 

. 168 





. 153 

. 136 

138 

. 131 

. 000 . 

185 

. 166 

. 152 

. 142 

141 

128 

. 189 

. 171 

156 

. 151 

. 145 . 

141 

. 000 

. 195 

. 173 

153 

154 


.137 .000 


94.9 

95 . 3 

95.5 95 . 

8 95.4 95.4 

94.9 

94 . 3 

95 - 5 95 - 7 95 . 8 

95 . 6 

95 . 0 

94.9 95 . 

4 95.8 95.6 

95 • 6 

95 . 6 

00.0 94.7 95.4 

95 .7 

95.7 

95.5 95 . 

3 94.6 95.4 

95 • 7 

95 . 7 

95 - 7 95.5 95•3 

95 . 8 

95 . 8 

95.6 95 . 

3 95.0 94.3 

95 . 5 

95 . 3 

9 A.fl 9 A . n 95 . A 

00 . 0 

95.6 

95.7 96 . 

0 95 - 8 95.7 




95.5 

00 . 0 

0 0•0 95 - 

8 96.0 95.8 

95.7 

95 . 6 

00.0 95.4 95.8 

96 . 0 

95.6 

95.5 00 . 

0 00.0 00.0 

0 0 . 0 

00 . 0 

0 0 . 0 0 0 . 0 0 0 . 0 

95.7 

96 . 0 

95.9 95 . 

9 95.4 95.1 

00 . 0 

95.8 


96 . 1 

96 . 1 

95.9 95 . 

6 95.5 0 0 . 0 

95.9 

96 . 0 

96.1 95.8 95.6 

00 . 0 

95 . 3 

96 -U 96 . 

1 95.8 95.7 

95.5 

00 . 0 


16.7 

10 . 1 

7 . 0 5.5 

4.8 4.1 3.8 

16.9 

10.2 

7.2 5.6 4.9 4.2 

17 . 0 

1 0 . 5 

7.4 5.9 

5 . 0 4.4 0 . 0 

17.3 

1 0 . 6 

7.6 6.0 5.1 4.6 

17.6 

1 0 . 7 

7.8 6 ■ £ 

5.1 4.7 4.5 

17.7 

11.3 

8.4 6.8 6 . 0 5.3 

18 . 0 

11.5 

S - 6 6 - 9 

6.2 5.4 0.0 

18.2 

11.8 

8.7 7.2 6.3 5.7 

0.0 1 

2.1 £ 

.9 7.4 t 

.4 5.9 0 . 0 1 

8.8 1 

3 . £ 

r * • 3 7-6 6 * • 7 6 . U 

0 . 0 C 

. 0 0 . 

0 0 . 0 0 . 

0 0 . 0 0 . 0 




19 . 1 

12.8 

9.9 3 . 1 

7.5 6.8 0 . 0 

19.5 

13 . 0 

10 . 1 8.4 7.7 6 . 


19.7 13.3 10.3 3.3 7.9 7.3 0.0 30.1 13.4 10.5 9.0 3.1 

8 577 4.91 

9 577 4.91 

10 577 4.91 

11 577 4.91 

12 577 4.91 

13 577 4.91 

14 577 4.91 

15 577 4.91 

16 577 4.91 

17 577 4.91 
13 577 4.91 

19 577 4.91 

20 577 4.91 
IS 2335 4.91 
13 2335 5. 

13 2335 5.5 
IS 2335 6. 

13 2335 6.5 

18 2385 6.37 
18 2885 7. 

18 2885 7.5 
IS 2835 8. 

18 2335 3.5 
18 8885 8.33 


1? 

1 154 

o _ 33 

18 

1 2 i'i fi 

8.83 

18 

14 0 0 

8.83 

18 

1 6 0 0 

8.33 

18 

1731 

8.83 

IS 

1 8 0 0 

3.83 

18 

£ 0 0 0 

8.83 

18 

2200 

3.83 

18 

2308 

8.83 

18 

24 00 

Q m 33 

18 

2600 

3.83 

18 

2800 

3.33 

13 

2885 

8.83 


♦EOR 

♦EDE 
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. 138 
. 176 
. 196 
. 125 
. 154 
. 0 00 


. 000 

. 147 



4. 0 
4.2 
4.9 
0. 0 
0. 0 

8 0 . 0 

.3 0. 









Output from CGEN 


PPM *k 

V L-N 

PMP CMM) SP 

•IT LE--KM 

EFF VOL 

KCU IN 

Sp’mT 

(FP3M PL 8 .'• =. 

4.910 

££65851858366 

. ££4 

94.900 

16.700 

9 . 

SP MT 

(FPOM PLG>=. 

4.910 

£ 099506545953 

.£09 

95.100 

13.400 

1 0 . 

SP MT 

577 

(FPOM PLG> =. 

4.91 0 

196 06647394c’ 

. 194 

95.300 

1 0 . 1 00 

11 . 

SP MT 

(FPOM PLG>=. 

4.910 

1S4£65£0 09665 

.177 

95.400 

8.550 

1 £. 

SP MT 

(FPOM PLG ' 1 = . 

4. 91 0 

174 084££74673 

. 160 

95.5 0 0 

7. 000 

13. 

SP MT 

(FPOM PLG> =. 

4. 91 0 

165191846 06 0£ 

. 150 

95.65 0 

6.850 

14. 

SP MT 

•FPOM PLG>=. 

4.910 

1573433888831 

. 140 

95.800 

5.500 

15. 

SP MT 

CFPOM PLG>=. 

4.91 0 

15 0354 0££5909 

. 137 

95.6 0 0 

5. 150 

IS. 

SP MT 

(FPOM PLG > =. 

4.91 0 

1440813001396 

. 134 

95.400 

4.800 

17. 

SP MT 

(FPOM PLG> =. 

4. 91 0 

1384134517809 

. 1 £9 

95.400 

4.45 0 

18. 

SP MT 

cr"7"?' 

(FPOM PLo'=. 

4.910 

1338614447906 

. 184 

95.400 

4. 100 

13. 

SP MT 

cr “7 -7 

(FPOM PLG> =. 

4. 910 

1£85533911983 

. 1 £3 

95.150 

3. 950 

30. 

SP MT 

(FPOM PLG>=. 

4. 91 0 

1£4830546401 6 

. 1 ££ 

94.900 

3.800 

18. 

SP MT 

(FPOM PLG> =. 

4.91 0 

£ 096 07£ 0 04393 

. 147 

95.500 

4.700 

18. 

SP MT 

(FPOM PLG."' — • 

5. 000 

£ 09131575 0777 

. 147 

95.500 

4. 760 

18. 

SP MT 

(FPOM PL 0 • = . 

5.50 0 

£06571£856301 

. 145 

95.500 

5.091 

18. 

SP MT 

£885. 

(FPOM PL 8 >=. 

6 . 000 

£041368995148 

. 144 

95.5 0 0 

5.4 £3 

18. 

SP MT 

£885. 

(FPOM PLG > =. 

6 . 5 0 0 

£ 018107499669 

. 143 

95.500 

5.755 

18. 

SP MT 

£885. 

(FPOM PLG'=. 

6 . 870 

£001586730145 

. 14£ 

95.500 

6 . 0 0 0 

18. 

SP MT 

£885. 

(FPOM PLG> =. 

7. 00 1.1 

199581££74473 

. 148 

95.500 

6 . 086 

18. 

SP MT 

(FPOM PL G '• =. 

7. 5 00 

1974390167889 

.140 

95.5 0 0 

6.418 

18. 

SP MT 

£885. 

(FPOM PLG > = . 

8 . 000 

1953735683996 

. 139 

95.5 0 0 

6.749 

18. 

SP MT 

£885. 

(FPOM PLG> =. 

8 . 500 

1933780834534 

. 138 

95.500 

7. 081 

18. 

SP MT 

£885. 

(FPOM PLG>=. 

8.830 

19£0963068465 

. 137 

95.500 

7.300 
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RP M * 

< V L-f 

1 ’WF (►•!■») S 

C WT LB/<M 

EFF 

V^L XC'I I'l 

18. 

1154. 

8.830 

. 131 

95.100 

8.800 

SP MT 

■: FPOM .9LG ■ 

■=.139S£05175885 




18. 

1 £ 0 0. 

8. S3 0 

. 131 

95.132 

8.8 0 0 

SP MT 

(FPOM 8LG) 

' =.141015 0£15183 




IS. 

14 0 0. 

8. 830 

. 130 

95.871 

6.800 

SP MT 

(FPOM ftLG 

' = .147078 08£0914 




1S. 

16 0 0. 

8. 830 

. 1 £9 

95.409 

6.800 

SP MT 

(FPOM 9LG 

•=.15314114S6845 




18. 

1731. 

8. 830 

. 188 

95.500 

6.800 

SP MT 

(FPOM ,9LG. 

'=.15711£4473398 




18. 

18 0 0. 

8. 830 

. 130 

> 95.484 

6 • 86 0 

SP MT 

'FPOM ,9 LG 

■ ~. 15 9 £ 0 4 £ 0 3 £ 3 7 5 




IS. 

8 0 0 0. 

8. 830 

. 134 

95.380 

7. 033 

SP MT 

(FPOM 8LG 

1 =. 18 5 £ 8 7 £8381 08 




18. 

££00. 

8. 830 

. 139 

95.£58 

7.806 

SP MT 

(FPOM 8LG. 

> = .17133 03843837 




IS . 

£308. 

8.83 0 

. 141 

95.£00 

7.300 

SP MT 

(FPOM 8LG 

>=.1746043770932 




18. 

£4 00. 
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Fig 11. Comparison of Algorithm and Subroutine, Changing RPM 
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